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Abstract 

In an earlier paper (J. Phys. A: Math. Gen. 26 (1993) 4689) we 
introduced the notion of cycle percolation in the Bose gas and conjec- 
tured that it occurs if and only if there is Bose-Einstein condensation. 
Here we give a complete proof of this statement for the perfect and the 
imperfect (mean-field) Bose gas and also show that in the condensate 
there is an infinite number of macroscopic cycles. 
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1 Introduction 

It was perhaps Feynman who first emphasised the importance of long per- 
mutation cycles in the description of the A-transition in liquid helium and, 
more generally, of Bose- Einstein condensation (Feynman 1953). Some years 
ago the present author gave a mathematically precise formulation of this idea 
by introducing the notion of cycle percolation in Bose systems (Siito 1993, 
hereafter referred to as I). In I it was suggested that this concept could serve 
as an alternative characterization of Bose-Einstein condensation (BEC). The 
examples of the ground state of a system of bosons and the three dimensional 
ideal Bose gas were treated. For the latter it was shown that BEC indeed 
implies cycle percolation, but the proof of the absence of percolation in the 
absence of BEC was missing. Since then, this kind of description has gained 
some new interest (e.g. Bund and Schakel 1999, Schakel 2000, Martin 2001, 
Ueltschi 2002). Therefore, it seems to be worthwhile to complete the earlier 
proof and to make the claim stronger by showing the occurrence of an in- 
finity of macroscopic cycles in the condensate. As it will become clear, our 
conclusions also apply to the imperfect Bose gas. 

Let us recall the principal definitions and relevant results of I. Let A be 
a cube of side L of the d dimensional Euclidean space, and denote Hi^^n the 
Hamiltonian of (interacting) bosons in A, taken, e.g., with periodic bound- 
ary conditions. By making explicit the summation over the permutations of 
A^ particles, necessary for symmetrization, the corresponding canonical par- 
tition function reads 



Here Sn is the group of permutations of {1, 2, ... , A^} and U{g) is the unitary 
representation of g in the A^-particle Hilbert space where Ti = L2(A). 
In the second line of ([T]) the summation goes over the conjugation classes 
of Sn- Each g in the same class giving the same contribution, see (1.2.5), 
g{nj} can be any permutation with rij cycles of length j. Because each term 
contributing to the sum is positive, cf equation (1.2.18), 
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can be interpreted as the probability of the occurrence of the permutation 
g in the canonical ensemble. Then, we can introduce random variables, for 
example ^i, the length of the cycle containing 1 (i.e. particle no 1) and rij, 
the number of cycles of length j. Constancy of the probability ([2]) within a 
conjugation class means that the probability is sensitive only to the number 
and length of the cycles and not to their contents. Thus, 1 is assigned to any 
cycle with a probability proportional to the length of the cycle, which yields 
the conditional probability 

^A,iv(6 = iWj = m) =mj/N (3) 

and, as a consequence, 

m 

Clearly, J2f=i P'A,Ni^i = j) = 1- However, if we first perform the thermody- 
namic limit L oo, N oo in such a way that N/L'^ — p > 0, and then 
sum over the limiting probabilities -Pp(^i = j), we can state only that 

oo 

E^p(^i=^')<l- (5) 

We speak about cycle percolation, in an obvious analogy with site or bond 
percolation, if in this relation the strict inequality holds. In I we conjectured 
that this happens if and only if there is BEC. Half of the conjecture, namely 
that BEC implies cycle percolation was shown for the free Bose gas. Now we 
give a complete proof of a somewhat stronger assertion. 

In the theorem below by condensate we mean the ensemble of particles 
in the one-particle ground state, when BEC occurs. Throughout the paper 
we work with periodic boundary conditions, so these are particles with zero 
wavevector. Furthermore, a macroscopic cycle is a cycle of positive den- 
sity, i.e. a cycle containing a non-vanishing fraction of the total number of 
particles. 

Theorem In the perfect and mean-field Bose gases Bose-Einstein conden- 
sation occurs if and only if there is cycle percolation. In the condensate the 
number of macroscopic cycles is (countable) infinite. Moreover, there are no 
finite cycles in the condensate and no macroscopic cycles outside it. 
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In the following section we prove the theorem in three steps. First, we 
recall the proof presented in I for the existence of cycle percolation when 
there is condensation. The second and third parts contain the new results: 
the absence of percolation in the absence of condensation and the analysis of 
the cycles in the condensate. 



2 Proof of the Theorem 



2.1 Proof of cycle percolation in case of condensation 

For the ideal Bose gas 

i=l 

the N-particle kinetic energy operator with, as we suppose in this paper, 
periodic boundary conditions in A. In rewriting ([1]), we make use of the fact 
that f/ is a representation of and, therefore, U{g{nj}) can be decomposed 
into a product of '^rij commuting factors, each corresponding to a cycle. 
Because now exp(— /^Ta^at) also factorizes, the trace will factorize according 
to the cycles, and the partition function reads 

with the trace over 7i. The key to the proof is provided by the bounds 

^ iV < Tre-^'^^-'^ < { — ^ + iV . (8) 



These are direct consequences of the inequalities 

applied with a = vrjA^/L^ where = h^j2ii[3/m, the thermal de Broglie 
wavelength. Now equation (jl]) with ([7j) yields 

PKAii = j) = ^^T^ Tr e-^^^'^- < i Tr e-^^^^^^ . (10) 
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The inequality holds because Qa^n-j < Qa,n for periodic boundary condi- 
tions. (See (1.6.16) for Dirichlet boundary condition.) Taking the thermody- 
namic limit and using the upper bound in ([8]), we find 



1 



(11) 



This leads to 
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where 



oo 
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a 



(13) 



According to the well-known condition (e.g. Huang 1987), there is BEC in 
the free Bose gas if and only if the right-hand side of (fT2l) is less than 1, 
which implies cycle percolation. 

2.2 Proof of the absense of percolation if there is no 
condensation 

Here we refer to the strong equivalence of ensembles in the absence of BEC. 
As we show below, for any fugacity z < 1 and any fixed j the grand-canonical 
probability distribution of nj/L'^ becomes degenerate in the thermodynamic 
limit, i.e. will be concentrated onto a single value. As a consequence, the 
canonical probability distribution of rij/L'^ will also be degenerate and con- 
centrated on the same value, provided that the particle density p corresponds 
to z. So it will be possible to compute Pp from the grand-canonical proba- 
bilities Pz- First we show the absence of cycle percolation in the absence of 
BEC in the grand-canonical ensemble. 
The grand-canonical partition function 



z^Qa,n =J27J}Y1 Tr^^f/;v(^)e 
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can also be considered as the generator of a probability distribution over 
USn. Since any g & is naturally embedded into S'at/ for any N' > N (by 
setting g{i) = i for i > N), it is more precise to say that now probabilities 
are assigned to the representations Un^q), 

PaAUn{9)) = j;^J'n-UN{9)e-^'''^^ (15) 

or, as the second line of (fT^ shows it, to terminating sequences of non- 
negative integers, each {uj} with ^ n j < oo corresponding to a conjugation 
classe of for = 'Y^juj. As in the case of the canonical ensemble, we 
will consider Uj to be random variables. For a general Hamiltonian any finite 
number of them has a coupled joint distribution. However, if for the ideal 
Bose gas we replace ([7]) into equation f|T^ . we arrive at 



(16) 

showing that the Uj are independent random variables with probabilities 

Pa,.(% = n) = exp |-y Tr e-^^'^'^'^ | ^ Tr e-^'^^^'^^ . (17) 

The L-dependence of (fT7|) leads to the asymptotic degeneracy of the grand- 
canonical distribution of Uj/L'^ for z < 1: 

limPA,.(7^<x) = e(x-p,) (18) 

where G is the Heaviside function and 

P. = -TT^ • (19) 

Equation ffTSjl is the manifestation of the following large deviation principle. 
Let < a < 6, then 

1 TL ■ 

lim — In Pa,, (a <-±<h) = - inf /,(x) . (20) 

L—>oo i^" h a<x<b 

where the rate function Ij is given by 

X 

Ij{x) = pj — X + xln — . (21) 

Pj 
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To prove (120]) and (1T51) . we extend the right-hand-side of (1T7|) to real positive 
values of n and define 

/.W=exp|-^TVe-'"-»}^(^Tle-'™..)" (22) 

for X > 0. Then, by using (IHD and Stirling's formula, asymptotically, as 
L — > oo, we find 

" TISP^""'" ■ '''' 

Equations f|T8|) and fl20l) follow from fl23|) by noting that Ij{x) is strictly 
convex for x > and Ij{pj) = is its minimum. For the distribution of nj this 
implies that asymptotically nj is concentrated to an 0(L'^/^)-neighbourhood 
of pjL'^. 

As in the canonical case, is the length of the cycle containing 1. In 
finite volumes 

Pa,z{^i = j) = ^ PaA^i = = mi,i> l)PA,z{ni = mi,i> 1) 
= E ^^^^,UP^,ML' = mjL'^). (24) 

With our large-deviation result then 

P-Mi = j) = lim Pa,.(6 = j) = ^ = , . . (25) 

Since the denominator is finite for z < 1, we find Xljli -Pz(Ci = j) = 1 indeed. 
Now for z < 1 p and z are related through 

gm{^) = pA| . (26) 



Substituting this expression into (]25ll we obtain Pp and the absence of cycle 
percolation in the canonical ensemble. 
Two remarks are in order here: 

1. As we have mentioned at the beginning of this section, the asymp- 
totic degeneracy f|T8l) implies that the canonical distribution of nj/L'^ is also 
asymptotically degenerate. Knowing this, we can derive Pp{^i = j) with- 
out passing through Pz{C,i = j). Indeed, in equation (jlj) for j = 1, . . . ,iV 
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the averages {nj)\^N can be approximated by the set {rij} belonging to the 
largest term of ([7j), which can be found by conditional maximization. The 
asymptotically valid result is 

^ J^^iWb ^^^^ 

Here appears via a Lagrange multiplier and satisfies 



N 

Z 



E 



ik = P^B. (28) 



If p < pciP) = 5frf/2(l)/A^ then z^ < 1- Writing equation for and 
and taking the difference one can see that for any fixed p > 0, z^ > ^at+i. 
So for p < pc zn I z < 1 which, hence, is just the solution for z of equation 
( |26|) . In sum, 

P^ACi = j) X - Pp(ei = j) = -TT^TT • (29) 



2. Comparing the expression fl2Sl) of Pz=i{C,i = j) with the right member 
of the inequality (ITT!) , we see that the inequality saturates at p = pc- In fact, 
for p > Pc there is equality in (ITT!) : Although for p > Pc there is no asymptotic 
equivalence between the canonical and grand canonical distributions of nj/L'^ 
(reflecting the fact that the canonical and grand canonical Gibbs states are 
related through the nondegenerate Kac density, see Cannon 1973 and Lewis 
and Pule 1974), the former still becomes asymptotically degenerate. So the 
strong equivalence of the canonical ensemble with a microcanonical one (with 
rij fixed) persists. The largest term of ([7]) is given by fl27j) . From fl28|) it 
follows that zn > 1 and converges to 1, therefore 

Pp(6=j) = -^ if P>Pc. (30) 



2.3 Existence of an infinity of macroscopic cycles 

While Yl'jLi Ppi^i = j) < 1 means that infinite cycles occur with positive 
probability, this inequality does not inform us about the asymptotic size and 
number of large cycles. Below we show that in the domain of condensation 
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macroscopic cycles are in abundance, their number is infinite, and each of 
them is associated with the condensate. 

Let, therefore, d > 3 and fix a p > pdP) = gd/2{^)/ ^b- Thus, po = 
p — Pc > 0. Let, moreover, = N{L,p) = [pL"^]. Writing the trace in the 
basis of the eigenstates of Ta^a?, the partition function reads 

'5a,- = 4E^"'^^"^^"^"EW- (31) 



K 

7d ^. _ /fe2r,2 



Here K = {h, . . . , fc^v), h e A* = (27r/L)Z^ ek = {h^k^/2m), Nk{K) is the 
number of occurrence of k in the sequence K and gK = {kg-i(^i))f^i. From 
fl3T|) we can infer the probabihty of a pair ((?, K), 



PaA9, K) = (32) 

and that of K, Pa,n{K) = Y.g Pa Ad, K). 

The asymptotic distribution of the random variables was studied in 
detail by Buffet and Pule (1983). They showed that the canonical distri- 
bution of Nk/L''- becomes asymptotically degenerate, namely it converges to 
the Dirac delta concentrated at po if ^ = and at if 7^ 0. From their 
analysis it follows that ^1 = 0{N) is possible only if fci = 0. Moreover, let 

Al,, = {Ke (AX I \No{K)/L'' - Pol < e, Nk{K)/L^ < e for ^ O} , (33) 

then 

lim PA,7v(^L,e) = 1 for any e > . (34) 

We prove equation (34) by showing that the probability of the complement 
of Al^^ goes to zero. 

PkA^I,.) < PAA^Uk = oy) + J2PAAALAkY) 

where 

ALAk = 0) = {K: \No{K)/L^ - po| < e} 

and for k ^ 

AlM = {K : Nk{K)/L'' < e}. 
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Now PA,N{^L,e{k = Oy) goes to zero as L tends to infinity. With = e ^^^^^ 



E 



< X 



:j^d Z_/rrt= \eL 



E 



{ng}:Y:ng=N LLq 



UqX7 



because in the numerator we have a part of the sum of the denominator, 
over those sets of occupation numbers for which Y^Uq = N, Uk = and 
no > \^L'^] (observe that Xq = 1). But 

f3e{k) = = 7rA|nVL2 

where n = (rii, . . . , n^), = Yli=i ^1 with Ui integers. Therefore 

< 



exp{-7reA|L'^-V} 

n€Z''\{0} 

oo \ 

X: exp{-vr6A|L^-¥} -1 



l + 2^exp{-7reA|L'='-V} 
1=1 



with L going to infinity if ci > 2. 

The crucial property of the probabihties (l32l) is that they are the same 
for all permutations leaving K invariant. Due to this uniform distribution 
and to (IMI) . the discussion of the ground state cycle percolation in I becomes 
relevant and we can derive the analogues of equations (1.5.16) and (1.5.18) 
(notice the different meaning of c? - the number of spin states - in I). Related 
classical results can also be found in Feller 1968. 

We fix e < po/2 and choose a K in Al,,. Then Nq{K) > Nk{K) for 
every k 0. We consider a j such that No{K) > j > maxkjLo Nk{K) . Thus, 
cycles of length j can occur only among particles with wavevector 0. The 
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conditional expectation value of rij, provided that K is given, is 

K)^ = 2^mP^,^K -m\K)- \{g^S^\gK = K}\ " 

(35) 

Here |{- ■ ■}! means the number of elements of {■■■}• Now the number of 
permutations leaving K invariant is W^N^^iK.y'.^ so we obtain 

knj)K = {nj)No{K) = ^ i^., ^^|{^ e SN,{K)\nj = m}\ = ^ . (36) 

Indeed, for any M > j the uniform distribution Phiig) = 1/M! over Sm gives 
rise to 

{nj)M = ^mPM{nj = m) = ^ J^'^l'^fi' ^ SmIuj = m}\ = - ; (37) 

in order to see this it suffices to recall that the probability of a conjugation 
class {rii} of Sm (Xlj^'^i ~ ^) 



PM=n^{-) . (38) 

i=l ^ ^ 

From fl36|) we obtain 

Pa,^(6 - J|i^) - ^^^^^^.0 - ^ . (39) 

Only particles with the same wave vector can be in the same cycle. Therefore, 
ki = implies < Nq{K) and 

PaA^i > J\K) = (l - 4,0 . (40) 

Fix now < X < po/ p. Equations flM|) and PUI) imply 

hm Pa,^ f ^ > = ^ - X . (41) 
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Indeed, by choosing e < xp, 
lim Pa^n (% > x] = 



1-^)^5. (42) 

Po±e/ P 

The upper and lower bounds holding for all e < xp, we can send e to zero 
and find (gTj). 

Complementary to this result is 



Ppi^i =j\ki = 0)= lim Pa,jv(6 = = 0) = for any fixed j (43) 

L— >oo 

whose analogue also existed in the ground-state cycle percolation, see (1.5.13). 
Equation (H3|) can be obtained with the following argument. If we restrict the 
probability distribution P\^]\r{g\K) to the permutations of the Nq{K) parti- 
cles with wavevector, the restricted distribution is still uniform. Therefore, 
by (1371) . equation (l36ll extends to cycles of an arbitrary length among the 
zero wavevectors. As a consequence, if G Ai^e is chosen so that ki = 0, 
equation fl39l) also extends to all j < Nq{K) and together with flMl) implies 
fH3|) . Thus, the full probability fl30|) comes from ki ^ 0. 

The number of infinite cycles of positive density is (countable) infinite. 
The number of those with density larger than x < ^ can be obtained from 
equations (IMIl and (15^ : 

No 

lim V (n,)A,^ = ln^ . (44) 

j=\xN] ^ 

This number can be arbitrarily large if x is sufficiently small. For m > 
In p — In pq the expected number of cycles of density between e"'-'""^^-' and 
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„m+l„ a'rfl r, 

ln^^-ln^ = lne = l. 

P P 

The intervals [e~^™'^^\ e~"^) are disjoint, their number is infinite and on av- 
erage there belongs one infinite cycle to each interval. 

Finally, we note that the above results apply without any further ado to 
the imperfect Bose gas (cf Huang 1987). The Hamiltonian of this latter reads 

Ha,n = T^,N + (45) 

with some a > 0, so probabilities and averages in the canonical ensemble 
coincide with those of the perfect Bose gas. 
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